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In this work we study the Nambu-Jona-Lasinio model in the SU (2) version with repulsive vector coupling and 
apply it to quark stellar matter. We discuss the influence of the vector interaction on the equation of state (EoS) 
and study quark stars that are composed of pure quark matter with two flavors. We show that, increasing the 
vector coupling, we obtain more massive stars with larger radii for the same central energy density. 
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The SU(2) version of the Nambu-Jona-Lasinio 
model (NJL) [1] has been applied to study quark 
matter and to investigate quark and hybrid stars 
[2131415] . In this work, we will include a repulsive 
interaction in the original NJL in order to inves- 
tigate its effect in the quark stellar matter EoS 
and the structure of pure quark stars [6I7I8I. 

1. THE NJL MODEL WITH VECTOR 
COUPLING 

In this work we consider a two-flavor NJL La- 
grangian of the form ^ : 

C = qii-f^df, - m)q + gs[{qqf + [qiib^qf] 



-gviqi'^qf 



(1) 



where g is a fermion field with Nf — 2 flavors 
and Nc colors. Apart from the bare mass m, 
the Lagrangian is chirally symmetric {SU{2)l x 
SU{2)ii). We have considered interaction terms 
in the scalar, pseudoscalar-isovector and vector- 
isoscalar channels. The gs and gv are the scalar 
and vector couplings, respectively, and are as- 
sumed to be constants with dimension MeV^^. 

The mean field thermodynamic potential den- 
sity at zero temperature and > 0, that are the 
conditions inside compact stars, with quark en- 
ergy Ep = y/p^ + M^, is calculated in the Hartree 



approximation and is given by [3]: 
n{0,^i;M,t^R)^-2NfNc 
+ (fiR - Ep) X e{^ji - Ep)} 

4=gv 



(M - m)2 



4gs 



(2) 



The thermodynamic potential density de- 
pends on two parameters, the dynamical fermion 
mass M and the renormalized quark chemical po- 
tential /Zfl, which are related to the scalar (ipip) 
and vector (V'^V') densities at the chemical poten- 
tial n by [QflO] : 



M = m-2gs{^pi^), 



(3) 
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fi -2gv{'ip'^'ip). (4) 

Eq. ([3]) is the well-known NJL gap equation for 
the dynamical mass M . In general, this equation 
has more than one solution, the physical being 
the one that minimizes f2. The vacuum contri- 
bution to the thermodynamical potential density 
diverges and has to be regularized. This regular- 
ization is carried on by introducing a cutoff A. 

Once we have the thermodynamic potential 
density, other thermodynamic quantities can be 
calculated in the standard way. In particular the 
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baryon number density, energy density and pres- 
sure are given respectively by: 

= (5) 

e = n + ^iNcP, (6) 

and 

P = (7) 

where the energy density and the pressure of the 
physical vacuum was set to zero. Then, for a 
given baryon chemical potential ^, these expres- 
sions need to be evaluated for the values of M and 
fiR which minimize the thermodynamic potential. 

1.1. Stability of Quark Matter 

The Nambu-Jona-Lasinio (NJL) model is one 
of the most popular models for studying the spon- 
taneous chiral symmetry breaking and its restora- 
tion at finite temperatures or densities [9110'. 
The model contains quarks with a current quark 
mass which acquire a constituent quark mass 
M through the quark interaction. In the mean 
field approximation, this interaction produces a 
quark mass, via a condensation process of quark- 
antiquark pairs that spontaneously breaks the 
chiral symmetry. A study of the quark matter 
stability was done in [3] , and we analyze here the 
pressure and energy per particle as a function of 
the baryonic density [11 . 

The character of the chiral phase transition 
crucially depends upon the model parameters, 
particularly the vector coupling constant. For 
gv/gs < 0.75 the phase transition is of first or- 
der. However, when gv/gs > 0.75 it becomes of 
second order since for different values of /i, only 
the M ^ minimum exists and it moves contin- 
uously to the M = minimum. In the Figure 
[l]we present the energy per particle for different 
values of the vector coupling. We observe that as 
the ratio gv / gs increases the quark matter gets 
less stable. 

Quark matter inside compact star is neutral 
with respect to electric charges and color. More- 
over, the matter is in /3 equilibrium, i.e., all pro- 
cess of beta decay of quarks involving electrons 
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Figure 1. Energy per baryon as a function of the 
baryonic density p/po for different values of the 
vector coupling. 

and antincutrinos and their inverse process are 
likely to occur at identical rates. Thus, we must 
impose both /3 equilibrium and charge neutrality. 
Throughout this paper we only consider the latest 
stage in the star evolution, when entropy reaches 
the maximum and neutrinos have already diffused 
out. The neutrino chemical potential is then set 
to zero. For /3 equilibrium matter we must add 
the contribution of leptons to the stellar energy 
density and pressure. The relation between the 
chemical potentials of light quarks is determined 
by beta decay and can be written as jl2j : 

p.d = P'U+Pe- (8) 

For charge neutrality we must impose: 

Pe = ^{2pu- Pd)- (9) 

Eqs. © and ^ need to be satisfied together 
with the self consistent relations for the quark 
effective mass during the numerical calculation. 
The EoS for the leptonic sector is also needed, 
and is considered as a relativist degenerate Fermi 
gas. 

The effect of the vector coupling in the quark 
matter EoS is shown in Figure [S] We can see that 
for a certain energy density an increase of gv / gs 
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increases the pressure. The equation of state EoS increase of gv /gs the stars have more mass and 

becomes stiffer as expected since we know that larger radii, 

the vector interaction is repulsive. 
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Figure 2. Pressure as a function of the energy 
density for different ratios of the vector and scalar 
coupling constants gv/gs- 
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Figure 3. The star mass is plotted as a function 
of the central energy density for different choices 
of the vector coupling. 



2. QUARK STARS PHENOMENOLOGY 

The quark matter EoS obtained in the previous 
section are use to solve the Tolman-Oppenheimer- 
Volkoff equations. These equations are integrated 
from the origin, where M{R = 0) = and e{R ~ 
0) = Ec, with £c the central energy density, up 
to a point where r = R and the pressure is zero. 
For each EoS we obtain a unique relation between 
the star mass and its central energy density [TH 
113114115] . Figure [3] presents the curve M/Mq x 
Sc for different ratios gv/gs as a function of the 
central energy density, normalized by the mass of 
the Sun M© = 1.9891 x W^kg. As the vector 
coupling increases, we obtain more massive stars 
for the same central energy density. The reason 
lies in the repulsive vector coupling, which makes 
stiffer the EoS. The internal pressure increases for 
a given central energy and supports a larger star 
mass. In the Figure 2] we show the M/Mq x R 
diagram for quark stars. We see that with the 
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Figure 4. Mass-radius diagram obtained for four 
values of the vector coupling. 



The masses and radii of the maximum mass 
stable configuration obtained with the inclusion 
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Table 1 



Quark star properties for the EoS described in 


the text. 


gv/gs 


(M/Mq) 


i?(km) 


ec(fm-4) 





1.29 


7.36 


13.69 


0.25 


1.51 


7.99 


11.61 


0.5 


1.67 


8.50 


11.01 


0.75 


1.76 


8.79 


9.56 


1.0 


1.86 


9.11 


8.69 



of the vector couphngs in the SU(2) version of 
the NJL model are larger than the correspond- 
ing NJL without the repulsive term. The radius 
is still quite small, below 10 km, smaller than 
the average neutron star value, 12-15 km but the 
mass is larger than the value 1.2 Mq which many 
neutron stars show (see Table [Ij. 

3. CONCLUSIONS 

Our study only concerns pure quark stars with 
two flavors. At star surface {P = 0) the en- 
ergy per particle is greater than the nucleon mass 
and therefore within this model quark matter will 
hadronize. Thus, these results for pure quark 
stars need to be seen with caution and we may 
say that it is more likely to expect in NJL model 
hybrid stars with a quark matter core with two 
and also three quark flavors. A work in the SU(3) 
version of the NJL model with vector interaction 
is underway. A superconducting quark phase re- 
cently proposed should also be considered [16117) . 
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